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We show that there are solitons with fractional fermion number in integrable 
N=2 supersymmetric models. We obtain the soliton S-matrix for the minimal, N=2 
supersymmetric theory perturbed in the least relevant chiral primary field, the $(1,3) 
superfield. The perturbed theory has a nice Landau- Ginzburg description with a 
Chebyshev polynomial superpotential. We show that the S-matrix is a tensor product 
of an associated ordinary ADE minimal model S-matrix with a supersymmetric part. We 
calculate the ground-state energy in these theories and in the analogous N=l case and 
SU(2) coset models. In all cases, the ultraviolet limit is in agreement with the conformal 
field theory. 
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1. Introduction 



Integrable models have the striking property that in a collision all momenta are 
conserved individually, and that the n-body S-matrix factorizes into a product of two- 
body ones. The enormous number of constraints this implies means that the exact S- 
matrix can often be conjectured |I]J|. The exact S-matrix encodes, nonperturbatively, 
physical information about our theory; it is just a matter of extracting the information. 
One way to extract physical information from the exact S-matrix of a 1+1 dimensional 
integrable theory is to do thermodynamics via the thermodynamic Bethe ansatz [fH|]. If 
the S-matrix is exact, the results of the TBA are non-perturbative. For massive, integrable 
theories obtained from perturbing a conformal theory by a relevant operator, for instance, 
we can compare the Casimir energy predicted by the TBA with the Casimir energy as 
calculated in the perturbed conformal theory. Thus an important consistency check on a 
conjectured scattering theory is to see if the Casimir energy goes over to the correct central 
charge in the UV limit (where the perturbation goes away). The next corrections to the 
Casimir energy can also be computed from the conformal theory and compared with that 
predicted from the scattering theory. In this way, the S-matrix contains a tremendous 
amount of information about the theory. Some excited-state energies can also be obtained 
in this manner |j|-[7|. 

Consider, for example, perturbing the p-th unitary minimal (N=0) conformal theory 
by the least relevant operator, the $(1,3) field: 
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S P ^S P + XJ $ (1 , 3) ^x. (1.1) 

Perturbative evidence suggests that this theory has an infinite number of conserved currents 
|| . It was explicitly shown that there is at least one kinematic conserved current (beyond 
the energy-momentum tensor), which insures that there must be completely elastic and 
factorizable scattering. For A > the perturbed theory (|1 . 1|) gives a RG trajectory flowing 
from the p-th minimal model in the UV limit to the (p — l)-th minimal model in the IR 
limit. Since there is a non-trivial IR limit, there must be massless particles in the spectrum; 
this situation is discussed in 0. For A < the perturbed theory ( |1.1|) is massive and the 
integrability requires factorized scattering. The exact scattering theory associated with 



this perturbed theory was conjectured in [|T0| - |T2| . The ground-state energy associated with 



this conjectured scattering theory was considered using the thermodynamic Bethe ansatz 
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13] and [12 1 . It was verified to the first few orders in the UV limit (A — > 0) that these 



results agree with those obtained from perturbing the conformal theory. 

In this paper we will propose a scattering theory for the N=2 supersymmetric analog 
of ( |l.lp , i.e. the N=2 supersymmetric minimal models perturbed in the least relevant 
operator, the $(1,3) superfield. We assume integrability for these theories — this has been 
shown to lowest order in perturbation theory, and it is argued that the supersymmetry 
prevents any corrections ||14|| . These perturbed theories have a nice Landau- Ginzburg 



description [15) with a Chebyshev polynomial as the perturbed superpotential. 

The Landau-Ginzburg analysis indicates that these N=2 theories are closely related 
to their N=0 analog scattering theories ( |1 . 1[ ) . In fact, we will show that the S-matrix 
for the N=2 theory is a direct product of the S-matrix for the N—0 analog ( |1 . 1|) with 
a basic N=2 supersymmetric S-matrix. This structure was conjectured in |16|], as the 
simplest S-matrix consistent with the quantum-group symmetry of the model. Taking 
the nonrenormalization of the superpotential to hold even nonperturbatively, an assertion 
which has been supported in many ways, we know the exact superpotential and the exact 
soliton spectrum, including the masses of the solitons [T7[] . Thus we show that the solitons 



conjectured in [|16| do exist. However, we have a very different interpretation of the 



supersymmetric part of the S-matrix — in |16[ the supersymmetry is postulated to act 
nonlocally, whereas here it acts locally. 

The striking feature of a Chebyshev superpotential is that it has a set of degenerate 
wells with equal-mass solitons interpolating between adjacent ones. The bosonic part is the 
type of N=0 Landau-Ginzburg potential expected for theory (|1.1|) fTH|]B. This provides a 
nice way of thinking about the Landau-Ginzburg potential for the iV=0 minimal model. It 
also indicates why the soliton spectrum of the N=2 theory has its tensor product structure: 
the bosonic N=0 solitons become a supermultiplet in the supersymmetric theory. (This 
tensor-product structure has also been conjectured for the N=l case [p!9| , |20[| .) We will show 
that this behavior generalizes to the N=2 D-series and to the exceptional invariants. In 
fact, all the N=2 theories we discuss have a close relation to an iV=0 model. The Landau- 
Ginzburg soliton picture has been discussed for many N=2 theories [ZT|, so making this 
correspondence can provide intuition concerning many non-supersymmetric models. 

One fascinating feature we find is that the solitons have fractional fermion number. 
This phenomena was first observed some time ago [22[ in a two-dimensional ^> 4 theory 



1 This degenerate- well potential can also be predicted by studying the associated lattice model 
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coupled to a Dirac fermion. In the Landau-Ginzburg picture, our models have the same 
structure as that treated in [22]. However, our models are exactly soluble and we obtain 
the exact S-matrices. The models discussed in this paper have solitons with charges of ±|. 
In our next paper [ 2"5f , we will find the S-matrix for the N=2 minimal models perturbed by 
the most-relevant operator, as well as for models with generalized Chebyshev potentials. 
These theories have more general (but still rational) fractional charges. In particular, the 
charges are multiples of l/(k+ 1) in the most-relevant perturbation of the k-th minimal 
model. 

We will provide checks on the N=2 scattering theories by verifying that the Casimir 
energy obtained from the TBA goes over to the correct central charges in the UV limit. 



We will also perform this analysis for the analog iV=l discrete series [|19l , |20j and for the 
perturbed SU{2) coset models[|0|, since these fit nicely into our general picture. We will 
also calculate tr(— 1) F in the N=2 theories as a further check on the TBA equations. 



2. N=2 soliton structure 

Consider a N=2 theory with an effective LG description characterized by some 
superpotential W(X) (X is a chiral superfield consisting of a complex boson and a Dirac 
fermion). Since the bosonic part of the potential is |W(X)| 2 , the vacua are the points 
in the complex X plane where dW = 0. The solitons Xij are the finite energy solutions 
to the equations of motion connecting the i-th and j-th vacua: X(a = — oo) = X^ % \ 
X(a = +oo) = X^ (as discussed in |17|], not all such kinks are to be regarded as 
fundamental solitons). In the soliton sector corresponding to a soliton Xij with mass 
m and rapidity 6>, the N—2 superalgebra is [ B3[ , [: 



-^2 -^2 



Q+ = Qi = Q+ = Q- = {Q+, Q-} = {Q-.Q+} = o 
{Q+,Q + } = 2AW {Q_,Q_} = 2(AW/)* , (2.1) 

{Q+,Q_}=2me e {Q+,Q-} = 2me~\ 

where AW = W(X^) — W(X^). It follows from the above algebra that there is 



a Bogomolny mass bound m > |AW| [24yi7| . The basic supermultiplet irreducible 



representation of ( |2.1| ) for solitons saturating this mass bound (the only type of solitons 
which we consider) is a doublet consisting of solitons u{9) and d{9) with: 

Q~\u{9)) = V2^e e/2 \d(8)) Q + \u{9)) = uV2^e- B/2 \d{9)) 

Q+\d(9)) = V2^e e / 2 \u(9)) Q~\d(9)) = w*v / 2^e" e / 2 |w(6')), ( ""' 
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where ui = AW/m. All other actions annihilate the states. 

Our theory has a conserved U(l) charge F corresponding to fermion number. The 
generators have fermion number ±1, whereas Q have fermion number =pl (this 
notation is less bizarre at the conformal point where left and right fermion number are 
separately conserved). In soliton sectors the fermion number operator F generally picks 
up an additive constant piece, leading to fractional fermion number [p2| , |2"5| . |2"6|j . In fact, it 



can be shown via adiabatic or index theorem techniques P5] . |2T)[ | that the fractional part of 
the fermion number in a soliton sector X^j of our theory is given by 

f = -^-(Im\nW"(X))\^. (2.3) 

In our soliton doublet, u(6) has fermion number e and d(6) has fermion number e — 1, 
explaining why they haven't been labeled as boson and fermion. We will often label the 
solitons by their fermion number. The phenomena of fractional fermion number in 1+1 
dimensions occurs in physical polymer systems where one extra bond (fermion) can be 
distorted into n solitons and, thus, each soliton carries fermion number 1/n f27fl . 

The supersymmetry is defined on multi-particle states in the usual manner. Since 
Q is fermionic, one picks up phases when Q is brought through a particle with fermion 
number. For example, bringing Q through a fermion results in a minus sign. Since we will 
have fractional charges, we must generalize this notion to 

Q ± \e 1 e 2 ) = \(Q ± e 1 )e 2 ) + e ±i ^\e 1 (Q ± e 2 )), (2.4) 



where the action on one soliton (Q e) is defined as in (|2.2| ). The charges Q act with 
the same phases as Q 1 ^. In notation analogous to that used for coproducts in a quantum 
group, (|2.4j) reads 



ACQ*) = Q* ® 1 + e ±i7rF ® Q*, (2.5) 
where F is the fermion-number operator. This similarity with the quantum-group action 



is not a coincidence, since N=2 supersymmetry is a special case of a quantum group |16[] . 

The fractional fermion number is crucial for obtaining the correct soliton content and 
S-matrices. In |T7[], for example, the fractional charge structure was ignored. It was thus 



mistakenly conjectured (in the context of minimal models perturbed in the most relevant 
operator) on the basis of CPT that each soliton supermultiplet should be a quadruplet — a 
tensor product of two copies of the basic soliton doublet discussed above. Taking proper 
account of the fractional fermion number, one sees that this doubling is unnecessary. 
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Corresponding to this doubling, the S-matrix obtained in [17] is the tensor product of the 
correct S-matrix with its complex conjugate. The S-matrices for this class of integrable 
N=2 theories and the thermodynamic calculations (which confirm the S-matrices) will be 
discussed in [23|. In this paper we focus on another (simpler) class of integrable N=2 
theories, N=2 minimal models perturbed in the least relevant chiral primary field. 



3. Chebyshev Superpotentials 

There is considerable evidence that the N=2 minimal models perturbed in the least 
relevant chiral primary field (the $(1,3) perturbation) are integrable jl4]]. Using the results 
of [^§1 , the effective superpotentials characterizing these perturbed theories are given by 
Chebyshev polynomials i: 

WWX = 2cos») = ^±^, (3.1) 

so that W(X) = X k+2 /(k + 2) — X k + . . .. For convenience, we have set the perturbing 
parameter to one (the powers of this parameter can be put back in by charge counting; 
e.g., W5 = X 5 /5 — PX 3 + P 2 X). These perturbed theories are intimately connected with 
SU(2)k- The chiral ring structure constants (in the natural basis ||28|| ), for example, are 
the fusion rules of SU(2) k Q. 

The vacua of our theory (|3.1| ) are the k + 1 solutions of dW(X) — sin(/c + 2)6>/ sin# = 0: 

X (r) = 2cos-^- r = l,...,k + l. (3.2) 

k + 2 K ' 

Thus all the vacua are on a line in the X-plane. This (k+ l)-well potential is characteristic 
of the iV=0 analog theories ( |1 . 1|) as well — the main difference here is that we have also 
the fermions. Using the methods of ||17| , it is easy to show that there exists a fundamental 
soliton connecting each of the adjacent critical points (|372"|). Thus the spectrum consists 
of k solitons X r ( r +-\\ for r = 1, . . . k and their k antisolitons X^ r+ iy. Any other possible 
soliton will break apart into two or more of these solitons. Taking our solitons to saturate 
the mass bound m > \AW\, using the value of the superpotential (|3.1| ) at the critical 
points 

W{XW) = (3.3) 



2 The fact that our integrable perturbation should be in one of the flat directions of [ |28| | is 
supported by |2S||. Our results can be viewed as additional confirmation. 
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we see that our fundamental solitons connecting adjacent vacua all have equal mass 
m = \AW\ = 4/(fc + 2). 

Each of these solitons is a supermultiplet: the doublet discussed in the last section. 



From (|2.3|) , we find that u r (r+i)(0) has fermion number 1/2 and soliton d r ( r +i)(0) has 



fermion number —1/2. The corresponding 2k antisolitons have opposite fermion numbers. 



For the k = 1 case this structure is just that discussed in [^2J : the Dirac equation has one 
zero mode in the presence of a soliton. (The supersymmetry requires that there be at least 
one; it turns out that there is only one.) When quantized, this Dirac zero mode results in 
a doublet of states, and the charge-conjugation symmetry requires that they have charge 
±1/2. 

Our soliton spectrum of 4k particles, consisting of a two-dimensional supermultiplet 
( |2.2|) for each soliton connecting adjacent vacua, is exactly the structure conjectured in 
||16|| for these models; although, as we will detail in sect. 4, the action of supersymmetry 
is different. We will find an S-matrix for these solitons in the next section, and in sect. 6 
we will show that the TBA calculation gives the expected central charges in the conformal 
limit. If one ignores this supermultiplet structure, one obtains exactly the soliton structure 
conjectured [10, ill (and verified to a large extent [fL3lQ for the iV=0 minimal models (|1.1|) . 



3.1. D-series, Eq, E 7 



The theories ( |3.1[) correspond to perturbing the Ak+i minimal models in the least 
relevant operator. We found that the vacua of the theory correspond to the nodes of the 
Afc+i Dynkin diagram and the fundamental solitons correspond to the lines in the Dynkin 
diagram connecting the nodes. There are more N=2 minimal models, corresponding to the 



D-series and the exceptional models Eg, E? and Eg ||15|| . The D-series should be integrable 
when perturbed by the least-relevant operator, because each model is an orbifold of an A- 
series model treated in the previous subsection, and the conserved charge discussed in [|14 



is invariant under orbifolding. One expects that the exceptional models perturbed by their 
least-relevant operator to be integrable as well. The effective Landau-Ginzburg potentials 
for these perturbed theories are again obtained using the results of |p8||. For the D-series, 



Eg and Ej these potentials were obtained explicitly in pi 



WD k (x 1 y) = — — - cos(/c — H xy 2 where 2cos</> = 2 — x 

K J. Zi 

.3 4 2 

ry 01 01 H* 

w ^y) = Y + -4- xy2 + Y-Y2 (3 ' 4) 

x 3 xy 3 2 Ax 2 xy x 1 
W Er (x,y) = — + — x y + — ^- + ^7 + 



3 3 9 9 81 4374 



where, again, we have scaled the perturbing parameter to one. Minimizing these potentials, 
one finds that the vacua of these theories correspond to the nodes of the relevant Dynkin 
diagrams. We conjecture that the fundamental solitons are the lines in the Dynkin 
diagrams connecting the nodes. Just as in the Ak+i case, W alternates in sign between 
connected nodes and, thus, all solitons have equal mass. The N=2 structure again requires 
each soliton to consist of a (u, d) doublet with charges ±1/2. 



4. The S-matrices 

4.1. X 3 /3 - PX (sine- Gordon at its N =2 point) 

The perturbed theory (|3.1|) with k = 1 corresponds to the sine-Gordon model, with 
coupling at the N=2 point: (3 2 = |87r in the conventional normalization. The N=2 
symmetry is a special case of the general quantum group symmetry of the sine-Gordon 
model at any coupling [j32| , |33|| . Even though the X s /3 — (3X model is the "same" as the 
sine-Gordon at this coupling, we find the action of supersymmetry in the two descriptions 
is very different. The difference arises from the fact that sine-Gordon and the X 3 /3 — f3X 
are different local projections of the same theory fl34|. This is analogous to the Ising 



model, where the spectrum can be either a free fermion or a strongly-interacting boson 
with S- matrix S = — 1. They can be mapped on to each other by the non-local Jordan- 
Wigner transformation, but they most definitely are not identical. The same type of 
behavior happens in our N=2 model. The standard sine-Gordon S-matrix at this point 
|lj] does not have an obvious supersymmetry: it must be realized nonlocally on the sine- 
Gordon solitons pTD|JT^1 . Even though the S-matrix in this section is formally the same as 



that of sine-Gordon at the appropriate coupling, ours does in fact allow a local action of 
supersymmetry. It is likely that there is an analog of the Jordan- Wigner transformation 
for this model, but we have not found one. 

Fermion-number conservation means that the S-matrix for the process | J(9i)K(02)) — > 
\L(02)M(0i)) (corresponding to a soliton X12 colliding with a soliton X21) is of the form 



(4.1) 





where, again, u and d have charges ±1/2 respectively (and their antiparticles have opposite 
charges). Demanding that the S- matrix commutes with the supersymmetry generators 
with the action (|2.2| ) and (|2.5| ) requires the S-matrix elements to be 



Z{9) cosh - 



b = Z(9)i sinh - ( 4 - 2 ) 
c = Z(9). 

Crossing, unitarity and the stipulation that there be no extra bound states fixes Z{9) to 
be 

7(9) TT r2 (-^i+^ r (24i+i + i)r(^+^-i) 



1 ,i f°° dt sin9 



exp(- 



coshf v 4 t cosh 2 ^ 

It is crucial in this derivation that AW/m is 1 for the X12 and — 1 for the X21. Tensoring the 
above S-matrix with itself gives the S-matrix obtained in ]17| where, as discussed in the last 
section, the basic soliton supermultiplet was unnecessarily doubled. Our thermodynamic 
calculations will confirm the above S-matrix. 

This S-matrix is the same as that of the sine-Gordon model at (3 2 = |87r if we identify 
u and d with the soliton of that model and d and u with the antisoliton; the fermion number 
in the local description becomes the topological charge in sine-Gordon. For purposes of the 
ground-state thermodynamics done in sect. 6, this distinction is irrelevant; both S-matrices 
have the same TBA system and hence the same Casimir energy. However, we emphasize 
that the interpretations of the action of supersymmetry are very different. 

4-2. N=2 minimal models with $^3 perturbation 

As discussed in sect. 3.1, the soliton structure resulting from the Chebyshev superpo- 
tential is that of the corresponding iV=0 minimal model with additional N=2 structure. 
It does not immediately follow that the S-matrices are a direct product of the N=0 soli- 
ton S-matrix and an N=2 part, because the representation of the supersymmetry algebra 
depends on AW/m, and this can depend on which solitons are being scattered. For 
this reason, the S-matrices for the N=2 minimal models perturbed by the most relevant 
operator (to be discussed in Pep cannot be a tensor product. However, with the $^3 
(least-relevant) perturbation, AW/m just alternates between ±1 (i.e. AW/m = (— l) r for 
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^r(r+i)) an d, thus, all of the two-dimensional supermultiplets obey the same supersym- 
metry algebra as the X 3 — (3X model discussed in the previous section. From the point 
of view of the supersymmetry, every scattering process for all Chebyshev superpotentials 
looks the same. Thus the S- matrix for (|3.1|) is a direct product: 



srw = sr u (^^M, (4.4) 

where Sj? =0 is the S-matrix found in [ |10| , |i~l~|| for the massive scattering theory coming from 
the iV=0 minimal model SU(2) 1 <g> SU(2) k /SU(2) k+1 perturbed as in ([□]), and the S- 
matrix Sj^S] 2 for the supersymmetric part is the one discussed in the previous subsection. 



This direct-product structure ( f4.4| ) was found in []T9|j20|1 for the N=l supersymmetric 



case and in |T(| for the N=2 case. Again, the S-matrix (|4.1|) allows a local action of 



supersymmetry, whereas the sine-Gordon S-matrix allows a nonlocal action. Both choices 
give the same thermodynamics. The nice thing about the LG analysis is that it gives the 
Chebyshev potential and hence the soliton structure directly. 

The D k series, E Q and E-j have the same tensor-product structure, and the N=2 part 
is exactly the same: 

sr\e) = sr\e)®s^{e) (4.5) 

where Sj^S] 2 is the basic N=2 S-matrix, obtained in this section, describing how the 
supermultiplet elements u and d scatter. The iV=0 part Sq =0 is easy to find by using the 
correspondence between a factorizable S-matrix and the Boltzmann weights for a lattice 
model. In models with solitons, each vacuum corresponds to a "height" in an RSOS lattice 
model; for example, the iV=0 A n model S-matrix is equal to the Boltzmann weights for 
the critical RSOS lattice model A n _i|T5]]. We conjecture that the D n and E S-matrices 



are the Boltzmann weights for the critical RSOS models associated with the D n and E 



Dynkin diagrams 36 . 



5. Thermodynamic Bethe Ansatz 

The TBA allows us to extract the Casimir energy for our theory on a circle of length R. 
The only required input are the soliton masses and their two-body elastic S-matrix. The 
basic idea is to compute the minimum free energy of a system of solitons at temperature 
T = 1/R living on a circle of length L — > oo. This is accomplished by finding the allowed 
energy levels on the circle and then filling the levels so as to minimize the free energy 
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F. The TBA analysis is thus in the grand canonical picture, with a chemical potential jj, 



P,p7| (general chemical potentials enter in the calculation of excited-state energies [|7|]). Our 
system can be viewed as a two-dimensional system on a torus with length L corresponding 
to the compactified space and length R = 1/T corresponding to compactified time. It thus 
follows that the ground-state energy E(R) of our system on a circle of periodicity R is 
given in terms of our minimum free energy by 

E(R) = RF min (R,L)/L as L -> oo, (5.1) 

when fx = 0. 

Consider iV particles with rapidities 9±, . . . , 9n on the circle. The allowed wavefunction 
must be invariant upon bringing every rapidity 9k around the circle, through the others 
and back to its starting point. In other words we must have 

where T(9k\9k+i, ■ ■ ■ , 9k-i) is the transfer matrix for bringing particles with rapidity 9k 
through the others. Our allowed wavefunctions must satisfy (|5.2| ) for every allowed rapidity. 
In other words, the allowed wavefunctions must be simultaneous eigenvectors of the transfer 
matrices for bringing each rapidity through the others and the eigenvalues must satisfy ( |5.2| ) 
(the Yang-Baxter relation ensures that the transfer matrices commute). This condition 
quantizes the allowed particle energies Ek = nrik cosh 9u on the circle. We thus need to 
know how to diagonalize the transfer matrix in order to even write down the TBA system 
of equations. It is for this reason that the TBA equations have been derived for only a few 
theories with non-diagonal S- matrices |l^,|35|]. In the next subsection we will review the 



warmup case where the S-matrix is diagonal. In the following subsection we will discuss 
some aspects of our case where the S-matrix is not diagonal. 

5.1. The TBA equations for diagonal S-matrices 

When the S-matrix is diagonal, i.e. the only scattering is of the type 0(6*1)6(6*2) — > 
S a b(9i — 6 , 2)6(6 , 2)a(6'i), the transfer matrix in (|5.2|) is diagonal for the iV particle 
wavefunction describing particles of rapidity 9i with the i-th particle having definite species 
ai. The single- valuedness condition ( |5.2| ) then becomes simply: 

\[S aka] {9 k -9 J ) = l, (5.3) 



e im k sinh 8 k L 
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an interacting- model generalization of the one particle relation = 2nnk/L. In the large 
L limit fl5"T3] ) gives the distribution P a (9) of allowed rapidity levels for particles of species 
a in terms of the distributions pb(O') of rapidity levels actually occupied by particles of 
species b: 

2nP a (9) = mLcosh(e) +J^f d6' p b {e')(j) ab {e - 0'), (5.4) 

b •* 

where 

MO) = -I • 

It is convenient to define quantities e a {9) by 



PaiO) X a e~^ 



P a {6) l + A B e-«-W 



(5.5) 



where the fugacity A a = e^ aR . This definition reflects the fact that all particle species in 
this paper (bosons and fermions) have S aa (0) = — 1 and, thus, at most one particle of a 
given species occupies a given momentum level I. 

The distributions p a (9) are now chosen so as to minimize the free energy. This yields 
the TBA integral equations for the e a (9) in the presence of a general chemical potential 

EM- 

-^(^-^)ln(l+A b e-^^). (5.6) 
Using ( |5.1|) it can then be shown that the ground-state energy of the system can be written 



in terms of the e a as: 



£ A (i?) = -V— [ d6cosh8\ii(l + \ a e- e ^). (5.7) 

~ 27T J 



Having obtained the system ( |5.6j ) appropriate for a given scattering theory, one can 
compute the Casimir energy E(R) by (numerically) solving this set of coupled integral 
equations. The UV limit (m — > 0) and the IR limit (m — ► oo) of the Casimir energy can 
be obtained analytically and are discussed in sects. 5.3 and 5.4 respectively. 



3 We define the S-matrix elements so that free fermions scattering off of one another have 
S=— 1: this corresponds to ordering the creation operators by a particle's position in space. This 
definition enables us to use the Yang-Baxter equation and do the thermodynamics ignoring the 
fact that there are particles with non-zero fermion number. 
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5.2. Obtaining TBA Equations for non- diagonal S -matrices 



We will describe the transfer matrices in ( |5.2|) (one matrix for each particle k) in terms 
of a single transfer matrix @. Denote the transfer matrix for bringing a particle of type 
a and rapidity 9 through N particles and ending up with a particle of type b by T a ^{6). 
We are here suppressing the dependence on the rapidities 9\, . . .9^ of the N particles as 
well as the q N labels Cj for the species of the iV particles before the scattering process and 
the q N labels di for the iV particles after the process (q is the number of different particle 
species). The components of T a b(9) can be written in terms of the S- matrix elements as 



(T ab {0))% = E ( d ~ e ^ S t k d (0~0 2 )... S d k f_ iCN (9 - 9 
allfc 



(5A 



Pictorially, this is much easier to digest; it is 

di d 2 



d 



N 



a 







k 2 k N - X 


b 











Cl 



C2 



CN 



where each intersection is an S-matrix element. Using the fact that all of our S-matrices 
satisfy S^f(0) = —S ac 5bdi scatterings at zero relative rapidity permute the colliding particles 
(with a sign), it is seen that the matrix tr o T(0) = ^2 a T aa (9) satisfies 



— tr o T(0 — 9 k ) — T{9 k \9 k+1 , ... ,9 



k-l, 



(5.9) 



the transfer matrices appearing in (|5.2|) . We can thus concentrate on diagonalizing tr o T(0) 
for general 9 and then set 9 to the different 9 k at the end. The Yang-Baxter relation ensures 
that the tr o T(0) commute for different 9 and thus can be simultaneously diagonalized for 
all 9 by a ^-independent set of eigenvectors; only the eigenvalues depend on 9. 

We do not need to find the eigenvalues of tr o T(0) explicitly in order to do the 
thermodynamics; rather, we derive constraint equations so that we can define eigenvalue 
densities just like the rapidity densities. This problem is equivalent to diagonalizing a 
transfer matrix in an integrable lattice model, where the S-matrix elements play the role 
of the lattice Boltzmann weights. There are a variety of techniques to do this ||38|| . One 



We are grateful to N. Reshetikhin for this presentation. 
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which is particularly elegant and which is applicable to our models is the Algebraic Bethe 
Ansatz 



We will discuss this in appendix A. In the next section we apply these 
constraint equations to a variety of models. The result is that we will be able to obtain a 
TBA system of equations of the form (|5.6| ) ( |5~77|) containing some extra fictitious, zero-mass, 



particle species to account for the eigenvalue contribution to (|5.2|). 



5. 3. Casimir Energy in the UV Limit 

The free energy ( |5.7| ) can be found explicitly in the m a R —>■ limit [40 ; |J]- For zero 
chemical potential, it is given by 



E{m a R -> 0) 



X, 



+ %a 



-CI 



Va 



1+Va 



(5.10) 



where C(x) is Rogers dilogarithm functional 

1 



£(*) = - 



dy 



lay ln(l - y) 



y 



X, 



exp(— e(0)), and y a = exp(— e(oo)). The ground-state energy is found by setting all 



A a = 1. It follows from (|5.6| ) that the constants x a are the solutions to the equations 

Y[(l + X b x b ) N - b , (5.11) 



where N a {, = ^- j d6(j) a b(9). (We restore the chemical potentials because we will use them 
in sect. 6.3) The constants y a in ( |5.10| ) are nonzero only for those species a' with m a > = 0, 
where they are the solutions to 



Va' = JJ(1 + XbVb' 



(5.12) 



where b' also runs only over massless species. The UV limit Casimir energy ( |5.10j ) is to be 
compared with that of the conformal theory 



E(R) = ^(h + h-±), 



(5.13) 



where (h, h) are the conformal dimensions of the operator which create this state. For 
unitary theories, the ground state has h = h = in ( |5.13|) . We can check a proposed 
scattering theory by comparing ( |5.10|) with (|5.13|) . Higher order in m a R corrections to 
( |5.10|) can be numerically computed from ( pT6|) and compared to perturbation theory about 
the UV limit conformal theory. 
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5.4- Casimir Energy in the IR limit 

The Casimir energy in the IR limit can be expanded in a series in e~ maR . From ( |5.6| ) 
it is easy to show that for m a R large, 

e a (0) nm a Rcosh9-J2 N ab'Ml + h>yb') (5.14) 

b' 

Plugging this into (|5.7|) , one finds 



E(mR^ oo) -> -V q a [ d9— cosh6> e - mai?cosh *, (5.15) 

a J 271 



where 



q a = l[(l + Xb'yb') Nab ' 



V 

This provides a useful check on the results: q a must be the number of particles of mass 
m a because in the R — > oo limit we have a very dilute gas of particles, and the first term 
in this expansion of the free energy is given by one-particle contributions only. In the case 
where there are solitons, q a is not required to be an integer: it is given by B 1 ^ , where B 
is the number of iV-particle configurations. 



6. TBA for Minimal Models Perturbed in the Least Relevant Operator 

In this section we obtain and analyze the TBA system of equations for the N=2 
scattering theories discussed in sect. 4. We will also consider the analogous N=l case and 
the case of all perturbed SU{2)i ® SU(2)k/ SU(2)k+i coset models, since this involves little 
additional work. From the form of the S-matrix (|4.4j) we see that the iV=0 case enters 
as an ingredient for our N=2 case (and also for the N=l case [ P^| , pi]| ). The RSOS S- 
matrix conjectured in [!(]] to describe the iV=0 minimal models SU (2)i®SU (2)k/SU (2)k+i 



perturbed in the least relevant operator is believed to lead to a TBA system (|5.6|) which 



can be described by the figure [13| 



1 2 3 fc-1 k 

o- -o- -o- -o- <s> 



Each node in the figure represents a particle species for the equations (|5.6[); the k — 1 open 



nodes have m a = while the node labeled ® is massive. The elements cf) a b in ( |5.6| ) are zero 
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unless the nodes associated with species a and b are connected by a line, in which case 
0ab(^) = (cosh(6>)) _1 . In other words, (f) a b(8) = l a b4>{Q) where l a b is the incidence matrix 
for the above diagram. This TBA system of coupled integral equations was analyzed in 
fl3|1 ; the system was proven to describe the RSOS S-matrix only for k = 1 and k = 2, 
but there is much evidence that the correspondence holds for the whole series. We shall 
assume that this is true in this section. 

It is straightforward to take the UV and IR limits in the manner of sects. 5.3 and 5.4. 
With this (j) a b, we have N a b = \l a b- The solution to the equations (|5.11|) is [^7j|T3 



sin 



7r(a+l) 
fc+3 



- I. 



sin 



(6.1) 



fc+3 



To solve the equations (|5.12|) , we notice that this case corresponds to the above Dynkin 
diagram with the /c-th node removed. Thus the solution is 



sin 



Va = 



r(a+l) 
fc+2 



1. 



sin 



(6.2) 



fc+2 



Plugging these into ( |5.10| ) , and using the amazing dilogarithm identity |PD| , |4~1 



sin 2 ^ \ 2k it 2 



2 7r(q+l) 



sin 



fc+3 



k + 3 6 



along with £(1 



x) 



= C(1)-C(x) and C(l 
E(mR -> 0) 



1, yields 



" (—) (1 
\6RJ K 



6 



(6.3) 



(6.4) 



(k + 2)(k + 3) J 

This, of course, is the correct UV limit. The first corrections to (|6.4|) also agree with the 
predictions from perturbed conformal field theory ]l3||f ■ The infrared limit gives a nice 
check: the number of iV-particle configurations for large iV for the /c-well soliton structure 
is given by q N , where q = 2 cos is the largest eigenvalue of the incidence matrix Z a &. 
Using the y a in (|5.15|) gives precisely this result. 



6.1. N=2X 3 -(3X 

In this section, we apply the discussion of the appendix to obtain the TBA equations 
for our basic A^=2 S-matrix Sj^Z] 2 (|4.2| ). The eigenvalues of the A^-particle transfer matrix 
are obtained from (A. 10) to be 

N n a N a 



A^;M)=^n^-^)ri cosh ( 



r=l 



-) | J i sinh( 

r=n+l 



(6.5) 
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where Z(9) is given by (O) and where we define 

N R-n, N 0-0 N 0-0 

^^n(" in M^))-\n^ in M^)+(-irn cosh (V 1)) - (6 - 6) 

r=\ j=l j=l 

The y r in (|6.5|) are the N distinct solutions to (A. 9), which here becomes 



for any choice of n = 1, . . . , N. Note that, given (|6.7| ), A has no poles in 9 and is bounded 
at infinity and, thus, does not depend on 9. 

Unitarity of the S-matrix means that the eigenvalues (|6.5|) have magnitude one, which 



is required for (|5.2|) . The corresponding solutions to condition (|6.7| ) are of the form 
y = zq + iit/2 and y = z-^ — in/2, with zq and z^ real. Defining the distributions of 
such solutions as Po(zo) and -Fq^o) respectively, the log of (|6.7| ) gives 

2KPi(zi) = I d9 M (6.8) 
J cosh(zi - 9) 

where i = 0, 0. 

Condition (|5.2| ) yields the relation 

mL cosh 9 + Im^- In A(0) = 2nP(9) (6.9) 
d9 

for the distribution of allowed rapidities. Thus the level density for the real particles 
depends on the eigenvalue distributions, just as (|5.4| ) depends on the rapidity distribution 
for all the particles. Each of the 2 N eigenvalues A(#, {y}) is determined by specifying 
whether a given y r contributes a cosh(i(6* — y r )) or an isinh(i(6' — y r )) to the product in 
( |6~5|) . (Remember that a given y r appears exactly once in this product.) We thus define the 
density P~(zi) of z% which contribute a cosh to the product (n contributions in total), and 
P l + (z l ) for the remaining contributions. We thus have the relation Pi(zi)=P i (zi)+P^~(zi). 
With these definitions, (|6.5|) becomes 

7ra4lnA(0)= I d9' P (9')Im^-\nZ(9 - 9') 
d9 J d9 

+ (6 ' 10) 



2 ^—f / cosh(^ — 9) 
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Now from (}4.3|), (|6.8|), and the integrals 

cos6>t f dy 1 9 



dt- 



4 J cosh 2 (^) J 2tv cosh(y — 9) coshy 7rsinh6> 

we obtain 



Relations (|6.9|) and (|6.10|) then become 



i=i 



The relations ( |6.8|) and ( |6.11| ) are constraints on the particle densities just like 
Thus we can proceed identically to the diagonal case. We minimize the free energy with 
respect to p and subject to the constraints ( |BT8| ) and ( |6.11|) . Defining 

m _ e" £(e) ^ + (*0 _ 

P(9) l + e- £ W Pi(zi) i + e - e i(*0 1 ' 

for i = 0, 0, it is seen that we obtain a TBA system system of the form ( |5.6|) and ( |5 . 7p 
consisting of three species, one massive and two massless, with 4> a b(9) = (cosh 9)~ 1 l a b, 
where l a b is the incidence matrix for 

1 

o — <s> — o 

The open nodes in the figure correspond to the massless pseudoparticles (the densities) 
and the solid node corresponds to the soliton (p density) with mass m = |AW|. 

To take the ultraviolet limit of this TBA system, we use our analysis of the non- 
supersymmetric case. Since l a b is the incidence matrix for As, we can read off the x a from 
(|Q|): xq = x-q = 2 and x\ = 3. The y a can be found by removing the center node, leaving 
two unconnected nodes and yo = y^ = 1. Using the dilogarithm identity ( |6.3| ) yields 

E (mR _> 0) - - (^) . (6.13) 

Thus our TBA system has the correct central charge c = 1 in the UV limit. The infrared 
limit also gives the correct result: we have 2 N different A^-particle configurations. 
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6.2. N =2 minimal model with $(1,3) perturbation 

To find the TBA system for the entire N=2 series, we must diagonalize the transfer 
matrix ( |5.8|) associated with the S-matrix ( |4.4| ). Since the S-matrix is a tensor product of 
an iV=0 S-matrix with a supersymmetric one, we can diagonalize the two parts separately.il 
Each diagonalization has already been done: the N—0 minimal-model part was discussed at 
the beginning of this section, while the supersymmetric part was done in the last subsection. 
Each piece contributes their associated pseudoparticles (the open nodes in the diagrams); 
the tensor-product structure means that the pseudoparticles from each piece do not couple 
to ones from the other. However, both types of pseudoparticle couple to the same massive 



particle, because both contribute to the single equation (|6.9|) . Therefore, the TBA system 
has a single massive particle. It follows that the TBA equations are described by the N—0 
figure for Sf? =0 joined at the massive node to the basic N=2 part: 




O 

We do not know of any significance to the fact that this is a Dynkin diagram. 

We can now compute the UV limit of the Casimir energy. The x a are given by J5T 

Xq = Xtt = k + 1 

, N2 (6-14) 
x a = (a + l) 2 -l 

To find the y a , we remove the massive node, leaving us with the Ak-i Dynkin diagram plus 
two individual nodes. Thus the y a for a = 1 . . . k are given by ( |B~2|) , while yo = Vq = 1. 
There exists another amazing dilogarithm identity (used in the D n section of [ |37| | ) 
which gives us 

E(mR^0) ^ - (—) (6.15) 

in agreement with the central charge for the N=2 minimal series. The IR limit is 
correct as well. It can also be shown from the TBA analysis that the bulk contribution to 
E(R) (terms proportional to e.g. R or RlnR in the small R expansion) vanishes. This is 
to be expected since our theories are supersymmetric. We will now calculate tr(— 1) F as 
another check on the supersymmetric nature of our theories. 



5 We emphasize that the model is not a product of two models, since that would yield a direct 
product of S-matrices. In our situation, a particle can be thought of as a bound-state of a particle 
in the fc-th minimal model with one in the N=2 X s — X model. 
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6.3. tr(-l) F 

To calculate tr(— 1) F , we use the fact that with chemical potentials /i a , the free energy 
is equal to — T\iiti(e^ aNa ^ T ). Thus in a theory with a diagonal S-matrix, one can calculate 
tr(— 1) F by setting the chemical potentials \i a = z7re a T, where e a is the fermion number 
of a particle of species a. Finding the appropriate chemical potentials in a non-diagonal 
scattering theory is more subtle, because the diagonalization of the transfer matrix means 
that we are not working in a basis of particle eigenstates. However, the thermodynamic 
expectation value of a symmetry operator is well-defined on each state in this new basis 
because these operators commute with the Hamiltonian and hence the transfer matrix. 
Therefore, we need to find the eigenvalue of the symmetry operator associated with each 
eigenvalue of the transfer matrix. With some methods of determining transfer matrix 
eigenvalues, it is not obvious how to do this. For example, in the inversion-relation method 



used in |T3|], one finds the constraint equations for the eigenvalues without constructing 
the eigenstate. However, using the Algebraic Bethe Ansatz, we have an explicit expression 
(A. 5) for the eigenstates, making this identification simple. 

In our series of models, the particles (u,d) have fermion number (§, — \). Our 
eigenvector ijj (A. 5) thus has fermion number n — (N/2). This is the fermion number 
associated with each eigenvalue ( |6.5| ). By definition, we have 

/ P++ f P+ = N-n 

■' (6.16) 
p = N. 



To calculate tr(— 1) F , we set (iqR = fi^R = —in, and fikR = The remaining chemical 
potentials for the pseudoparticles l,...,(/c — 1) remain zero; these particles arise from 
diagonalizing the part of the S-matrix independent of fermion number. 

To finish the calculation, we need to calculate the x a and y a from ( |5.11|) and ( |5.12|) . 



The result is that with these chemical potentials, x a = y a . Specifically, Xk = yt = 0, 
= Vo = x q = Vq = 1) and the remaining ones are as in ( |6.2| ). The fact that x a = y a 
means that E\(R) = for all R — not only in the UV limit. Thus 

tx{-l) F = e- E{ - R)L = + 0(1). (6.17) 

Since we are working in the large L limit, we have proven only that the leading term, 
proportional to L, vanishes. This, of course, is required in a supersymmetric theory — 
( |6.17| ) gives the number of bosonic ground states minus the number of fermionic ones, 
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because the others pair up and cancel. (For the k theory, tr(— 1) F = k + 1.) Thus this 
calculation gives no new information, but does provide a check on our TBA system, as 
well as giving as simple example of the application of the excited-state TBA analysis 
to models with pseudoparticles. 



6.4- N=l minimal models with $(1,3) superfield perturbation 

The S-matrix for the scattering theory coming from the N=l minimal model 
SU(2) 2 ® SU(2)k/ SU(2)k +2 perturbed in the least relevant, supersymmetry preserving 
operator (the $(1,3) superfield) was conjectured in JT5| , |2"m to be 



sr =i (#) = sr u (o) ® Sfc=T(0), (6.18) 

where S^ : T 1 1 =S^s^ ) , the S-matrix for the perturbed tricritical Ising model. This was 
conjectured because it is the simplest S-matrix consistent with the appropriate quantum- 
group symmetry and the supersymmetry [2D]. The first part of this tensor product 



arises because the Landau- Ginzburg picture [|T5[ predicts solitons interpolating between 
the minima of a Chebyshev superpotential, just as in the N=2 case. The second part 
gives the model its supersymmetric structure, which acts nonlocally on solitons Jlpj. The 
origin of this soliton structure is unclear. Presumably there is a bosonization of this series 
analogous to the sine-Gordon description of the N=2 X s — X model, and the solitons 
interpolate between vacua of this field. 

To find the TBA system, we must diagonalize the transfer matrix (|5.8|) associated 
with the S-matrix (|6.18| ). As in the N=2 series, the S-matrix is a tensor product, and 
we diagonalize the two parts separately. Each piece is an iV=0 minimal model, so this 
diagonalization has already been done [H| , and is discussed at the beginning of this section. 



As in the N=2 series, the resulting two sets of pseudoparticles are coupled through the 
massive node: 

1 2 k- 1 k k + 1 

o- -o- -o- -®- -o 



It is simple to take the UV limit of this system using the analysis of the N=0 minimal 
models. Since we now have k + 1 nodes instead of k, the x a are given by (|6.1| ) with the 
substitution k — ► (k + 1). When the massive node is removed, we obtain a diagram of 
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k — 1 nodes along with a single disconnected node. Thus the y a for a = 1 . . . k are given 
by (|6.2|) , while yk+i = 1- Using the dilogarithm identity ( |6.3|) again, we find that 



/ 7T W3 12 



B(roH ^ 0) ^_J ^_ - ■__ ), (6 . 19 ) 

in agreement with the conformal results. 

We note that the k = 2 model in this series corresponds to the first model in the 
N=2 series. However, the S-matrix ( |K.18| ) for k = 2 is not identical to the S-matrix 
Q41|) . In fact, the particle content is different: the particle content here is not the soliton 
doublet discussed above but instead is that of the tensor product of two triple- well soliton 
potentials. This is equivalent to making the vacua the nodes and the solitons the links 
of the Dynkin diagram for As or D<±, depending on how one does the tensor product. 
However, one can recover the S-matrix elements (|4.1[) from ( |6.18[) by taking an orbifold, 
where one mods out by a symmetry. The orbifold procedure is a generalization of the 
Kramers- Wannier high- to low-temperature duality. This redefines the states as linear 
combinations of the previous states, hence changing the S-matrix. The new S-matrix can 
be found in the manner of lattice-model orbifolds [4^] , and one recovers the S-matrix Q4.1|) . 



This was done in to relate an iV=0 D-series S-matrix (the tricritical three-state Potts 



model) to one in the A-series. In fact, this can be done for every other model in the D- 



series, as in the lattice models of |43[]. The orbifolding should not change the ground-state 



energy, but will affect excited-state energies [44 



6.5. SU(2)i®SU{2) k /SU{2) k+l 

The general 577(2) coset case is the obvious generalization. The S-matrix here is |20 



s?r° (*) = m ® s?=° m , (6.20) 

It follows that the TBA system is described by the diagram 

1 2 fc-1 k k+1 k+l-1 

O O O <8> O O 

The x a are given by ( |B~ll ) with the substitution k — > k + I — 1. The y a for a = 1, . . . , k are 
given by ( |6.2| ) , while for a = k + l,...,k + l — 1 they are given by ( |6.2| ) with a —>■ a — k and 
k — > I. The result for the energy is in agreement with the central charge of the conformal 
theory c = 3(^ + ^-4±^). 
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7. Zamolodchikov's c-function and more 



The ground state energy E(R) described in the previous section in terms of the 
solutions to our set of coupled integral equations can be thought of as a Zamolodchikov 
'c-function' |45[] for our massive theory: E(R) = —7rc(mR)/6R. It is not obvious from 
the TBA equations that this 'c-function' monotonically decreases in the infra-red, though 
this is observed to be the case in the numerical and analytic analysis of all physical TBA 
systems checkedi. 



In [p9| another 'c-function' was obtained for the Chebyshev theories using the N—2 

is 



structure. It is interesting to compare these results. The c-function obtained in [£9 
quite simple: 

Q 7 

c k (mR = z) = -z — u(z, r = 2k/(k + 2)) (7.1) 

Zi CliZ 

where u(z, r) is a solution of the differential equation u zz + z~ 1 u z = sinh u (which can be 
transformed to Painleve III). The boundary conditions of u are specified by (for r < 2) 

u(z, r) ~ rlogz + const. + . . . as z — > 0. (7-2) 

The values for r in the solutions (|7.1[ ) are obtained from regularity conditions. The only 
k dependence is in the boundary conditions of u. It is easily seen from ( |7.2| ) that ( |7.1| ) 
yields the correct central charges in the UV limit z — > 0. It would be interesting to 
directly calculate the next UV corrections to the 'c-function' from perturbation theory 
around the conformal theory and compare with that predicted by the TBA approach and 
that predicted by ( [7.1|) . Comparing the IR limits (z — > oo), it can be seen that the two 
approaches give different results: the c-function ( f7.1|) has contributions corresponding to 
only odd numbers of particles, whereas the TBA approach includes contributions from 
any number of particles. It remains to be seen if one can relate these approaches by some 
projection in the manner of 0. If this were possible it would be quite interesting: we 
could work with an differential equation (and a well-studied one at that) rather than the 
TBA system of coupled integral equations (|5.6|). 



T.R. Klassen and E. Melzer, private communication 
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8. Conclusions 



We have found the exact S-matrix for the N=2 discrete series perturbed by the 
least-relevant supersymmetry-preserving operator. The nice thing about the N=2 non- 
renormalization theorems is that they leave little guesswork; here one does not have to 
conjecture the particle content. Since our results fit in so nicely with the iV=0 discrete 
series, we believe that this lends more support to the conjectured iV=0 soliton picture, and 
that much intuition can be gained by using the N=2 models to study the iV=0 models. 

Another interesting aspect is the appearance of solitons with fractional fermion 
number: ±| here, and any rational number in our next paper. There are only a few other 
integrable models with such structure, and those have only fractional charge ±^ |46 



Therefore, this analysis may be valuable for studying experimentally-realizable systems 
of fractional fermion number. It is probably not possible to measure the S-matrix or 
the Casimir energy directly, but one knows many other things about these N=2 systems 
which may be observable. For example, all the critical exponents of these systems are 
known exactly. 



Acknowledgements 

We are very grateful to N. Reshetikhin for many helpful conversations on the TBA, 
and for providing the presentation of sect. 5.2. We would also like to thank C. Vafa for 
valuable discussions, and T. Klassen and E. Melzer for their comments on the manuscript. 
K.I. was supported by NSF grant PHY-87-14654 and an NSF graduate fellowship, and 
P.F. was supported by DOE grant DEAC02-89ER-40509. 



Appendix A. Algebraic Bethe Ansatz for the 6-vertex Model 

Consider a scattering theory consisting of two species u and d with the two particle 
S-matrix: 

du ud uu dd 




(A.l) 



Our S-matrix elements are taken to satisfy the condition 



a{6)d{e) + b(e)l(e) - c{e)c(e) = o. 
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(A.2) 



Drawing a scattering vertex with right or up pointing arrows to denote species u and left 
or down arrows to denote species d, the scattering theory (|A.l] ) corresponds to the 6-vertex 
model. Condition ( |A.2|) is known as the free fermion condition in the six-vertex model|3S 



We will discuss the algebraic Bethe ansatz technique for finding the eigenvalues of the 
transfer matrix for theory (|A.1| ). Note that our basic N=2 scattering theory is of this 



type; we are here being more general with future applications in mind. 
The matrix T a b{9) discussed in sect. 4.2 will be denoted by 

Iab[U) -\C(6) D(6)J> [A - 6) 

C{9), for example, is the transfer matrix for bringing a particle of type d and rapidity 9 
through N particles with rapidities 9i and ending up with a particle of type u coming out. 
As discussed in sect. 4.2, we will be interested in finding the eigenvalues of the 2 N x 2 N 
matrix tr a T = A(9) + D(6). 



Define O = | YliLi d(@i)}- Notice from ( |A.1| ) that for any rapidities we have C{9)Vt = 0. 



The state O is also an eigenstate of the matrices A and D: 

N N 

A(9)0, = Y[b(e-di)fl, D(9)Q = Y[d(9-9i)n. (A.4) 

i=l i=l 

The algebraic Bethe ansatz is that all eigenvectors of the matrix A + D can be written as 

n 



r=l 



for some choices of the y r and n (it follows from (|A.1| ) that the vector ( |A.5| ) has 



n u 



particles and N — n d particles). As a result of the Yang-Baxter relation, T a b{9) satisfies 
the relation 

Si c '(9 - y r )T alb (9)T c , d (y r ) = T cd/ (y r )T a¥ (9)St! } d/ (9 - y r ) (A.6) 
This yields the following useful commutation relations for the transfer matrices: 

A(9)B(y r ) = <1 ^ t \B{y r )A{e) - ^4^)^) 
b[9 - y r ) b{9 - y r ) 



D{9)B{y r ) = - a ^-^\ B {y r )A{9) + C j^-^B(9)D(y r ) 
b{9 - y r ) b{9 - y r ) 
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(A.7) 



provided b(6 — y r ) 7^ 0. We have used the condition ( |A.2|) to simplify the equations. It 
follows from these relations and ( |A.4| ) that A{9) + D(9) act on the states (|A.5| ) as 



N N 



[A(9) + D{9)]^ = TT a( f Vr) (f[ b(9 - Bi) + (-l) n TT a(9 - + bad terms, (Ai 



where the 'bad terms' are vectors related to (|A.5|) by some B(y r ) — > B(9 P ). All such bad 
terms are proportional to [A(y r ) + (— l) n D(y r )]Q and thus vanish if 

fJl a(Vr - Bi) 

this determines the y r . Actually, since our S- matrices all have 6(0) = 0, the argument 
(using ( |A.7| )) that all bad terms are proportional to ( |A.9| ) breaks down if some of the y r 
are the same. We thus require the solutions of (|A.9| ) to have no repeating values of y r . 
Our desired eigenvalues of tr T(9) are thus given by 

n (0 — ) N N 

X ^ y) = U a b(e _ y ) (II W - °i) + (-!) n II ~< d - (a.io) 

r=l ^ ^ r > t=l i=l 

for any choice of distinct solutions y r and n of ( |A.9| ). 
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